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Abstract 



We investigate unambiguitics in the extraction of the d — u asym- 
metry from semi- inclusive production of pions in DIS. The role of 
p^ , several effects beyond the quark-parton model (QPM) which lead to 

^D I Np ^ N^ and N^ ^ N^ and may therefore cloud such an extrac- 

ts , tion is studied. The results are discussed in the context of the recent 

pi i' HERMES data. We find that the interaction of the resolved photon 

JL . with the nucleon significantly modifies the observed d — u asymmetry, 

flj ' The exclusive elastic production of p mesons plays a similar role for 

r^ , the large-z data sample. Our estimate shows a rather small effect of 

the spectator mechanism. Nuclear effects in the deuteron also look 
potentially important but are difficult to estimate quantitatively. 
^% ' Throughout the paper we make in addition several general remarks 

!^ , concerning semi-inclusive and exclusive production of mesons. 



1 Introduction 

Since the NMC publication Q on the Gottfried Sum Rule violation the ef- 
fect oi d — u asymmetry was in the 90's one of the most intensively discussed 
problems of nucleon structure. The effect is clearly of nonperturbative na- 
ture and was qualitatively explained as due to the pion (meson) cloud in the 
nucleon (for recent reviews see Q). In order to shed more light on the na- 
ture of the Gottfried Sum Rule violation two different Drell-Yan experiments 



were proposed and performed |^, ^]. They measured the ratio cy^d /^p^ • 
The integrated result for the asymmetry from the more complete Fermilab 
experiment |^] is Jq [d — u] dx = 0.09 ± 0.02, to be compared with the 
NMC result: Jl[d - u] dx = 0.148 ± 0.039. The NMC asymmetry apears 
slightly bigger. It was suggested recently by two of us that the difference 
can be partly due to large higher-twist effects for the nonsinglet quantity 

It was proposed in Ref . Q how to use semi- inclusive production of pions 
to extract the asymmetry of light antiquarks in the nucleon. This method 
was applied recently by the HERMES collaboration at HERA [^. 

Recent results for semi-inclusive production of pions in polarized photo- 
production obtained at SLAC [^ have shown that spin asymmetry almost 
cancels for small transverse momenta of the outgoing pions, which seems to 
be another nonperturbative effect. This result was interpreted as a large 
VDM contribution Q for small transverse momenta. Only at large trans- 
verse momenta may the perturbative QCD processes reveal themselves and 
only then can one try to disentangle the polarized quark distributions in the 
nucleon. At low photon virtuality, as in the HERMES experiment where 
<Q^> ~ 2.3 GeV^, similar nonperturbative effects can be expected in the 
unpolarized case. 

In the present paper we try to examine the semi-inclusive production of 
pions in DIS as a source for measuring the d — u asymmetry. We investigate 
several effects, mostly of nonperturbative nature, which may modify the 
resulting asymmetry. In particular, making quantitative estimations, we 
focus on conclusions relevant for the HERMES experiment. 

2 Quark-parton model approach 

2.1 Extraction of the d — u asymmetry 

The most general five-fold cross section for one-particle semi-inclusive un- 
polarized lepton-hadron scattering can be expressed in terms of four inde- 



pendent semi-inclusive structure functions (see for instance [10|). If the 
azimuthal correlation between the lepton scattering plane and the hadron 
production plane is not studied, the number of independent structure func- 





Figure 1: Partonic mechanisms of pion production: (a) direct fragmenta- 
tion of quark into a pion (b) fragmentation of quark into an intermediate 
hadronic resonance and its subsequent decay. 



tions reduces to two. Then the cross section can be written as 

^™ [y''2xJ^i{x,Q'',z,pl,) + 2il-y)J^2ix,Q^z,pl,)] 



dxdQ'^dzdp'l j_ Q^x 



il) 



where x, y and Q^ are standard DIS variables, Ph,± is the transverse momen- 
tum of the detected hadron with respect to the momentum of the virtual 
photon and 



P-Ph 
P-q 



Eh 



(2) 



is a relativistically invariant variable which in the target rest frame is the 
fraction of the virtual photon energy v carried by the hadron. In the formula 
above, P, Ph and q are the four-momenta of the target nucleon, final hadron 
and virtual photon, respectively. 

If one is not interested in the transverse momentum distribution of the 
emitted hadron then the triple-differential cross section can be written in a 
more compact way 



da 



ina^ 



dxdQ'^dz Q'^x 



[y^2xJ^i{x, Q2, z) + 2(1 - y)J^2ix, Q^ z)] . (3) 



In the quark-parton model (QPM) only mechanisms shown in Fig.ffl are 
assumed. Usually in calculations one does not distinguish diagrams (a) and 
(b). It is commonly believed that diagram (b) can be included effectively on 



the same footing as diagram (a). We shall discuss later possible restrictions 
of such an approach. 

In the naive QPM the generalized semi-inclusive structure functions J-\ 
and T2 are related by the Callan-Gross relation leaving only one independent 
structure function, which can be written as 

T^^^{x, Q\ z) = Y, ejxqfix, Q^) ■ Df^^{z) , (4) 

/ 

where the sum runs over the quark/antiquark flavours / = u,d, s; qj are 
quark distribution functions and Df^T^(z) are so-called fragmentation func- 
tions Q. 

Quite a number of fragmentation functions can be reduced by the re- 
quirement of isospin symmetry and charge conjugation: 

DZ^{z)=Df{z)=Dl{z) = Dl{z) 

^D+{z) (5) 

for the favoured fragmentation and 

Df{z)=Df{z)=Dl-{z) = Df{z) 

^D^{z) (6) 

for the unfavoured fragmentation. For the strange/antistrange fragmenta- 
tion, in principle, a third type of fragmentation functions has to be assumed. 
In the following we shall assume simply: 

Df (z) = of (z) = Df (z) = Df (z) = Sred -D.iz), (7) 

where Sred is a reduction factor with respect to nonstrange quarks/antiquarks. 
In calculations we shall use Sred = l-Of] 

Now, in the quark-parton model (^ using symmetry relations (^g) one 
can combine semi-inclusive cross sections for the production of positive and 
negative pions on proton and neutron targets and isolate a quantity sensitive 
to the flavour asymmetry Q 

d{x) - u{x) _ J{z)[l - r{x, z)] - [1 + r{x, z)] 

u{x) - d{x) ~ J{z)[l - r{x, z)] + [1 + r{x, z)] ' ^ ^ 



^This appears to be not very important in practice. 



where J{z) = l l+^-ii^^+li and r{x,z) = ^^ (-■^)-^^^ (-■-) ig ^ ratio of 

differences of charged pion yields on proton and neutron. It is straightfor- 
ward to see that the fragmentation of strange quarks/antiquarks cancels in 
the quantity r(x, z). It is also worth noting that the r.h.s. in Eq.(P), formally 
dependent on z, gives a quantity independent of z. 

Thus semi-inclusive production of charged pions in DIS allows us to 
determine the asymmetry of light sea quarks. Although this is on condition 
that the QPM works well, that is, one may neglect the influence of other 
possible mechanisms. 

2.2 Intermediate resonances in the fragmentation 

It is a well known fact that pions produced directly in the fragmentation 
process constitute only a fraction of all pions registered in detectors. The 
contribution of pions coming from the decay of heavier mesons is of the same 



order of magnitude [|1^, 13 , |T4| . 

Because the intermediate resonances originate from the fragmentation of 
the struck quark their contribution can be included into an overall "effec- 
tive" fragmentation function. Modern analyses of fragmentation functions 
do not treat intermediate resonances explicitly, just include them effectively 
by fitting total inclusive data. However, it is not clear, a priori, whether 
under a more detailed consideration such an effective treatment is correct 
and whether resonances do not disturb the identity (|g). 

To treat the intermediate resonances explicitly we write down the frag- 
mentation function as a sum of two parts: a direct fragmentation component 
(Figjl|a) and a resonance component (Fig.||b): 

Df^^ = Df^^ + Y, Df^R^^ , (9) 

R 

where Df^T^ is a fragmentation function of the direct fragmentation of a 
quark / into a pion vr, Df^R^T^ describes the production of a pion vr through 
an intermediate resonance R and the sum runs over all possible resonances. 
It is known experimentally that for pion production the vector meson inter- 
mediate states are the most important. For not too small z, neglecting for 
simplicity transverse momenta, the contribution of the resonance R to the 
fragmentation function can be approximated as 

Df^n^^iz) = f Df^Riz') ■ Ir^^ (^) dz , (10) 



where zo = max{z, z^^^), z^^^^ is the minimal possible z of the resonance R, 
Df^fi is a fragmentation function for the direct fragmentation of a quark 
/ into a resonance R, and fn^n describes the decay of the resonance R to 
pionic channels. 

The fragmentation process transforms quarks with the third component 

of isospin /| = ibi <^ s^-^ into measured pions with /? = ±1 <^ 



_ > i.e. there are two initial and two final states of fragmentation with 

TT J 

respect to I^. If the quark hadronization is driven by the strong interaction 
(isospin symmetric) then, in the case of direct fragmentation, one naturally 
obtains only two kinds of fragmentation functions related by (0-0). For 
the resonance contribution D^^^ = YliR^f^R^-^^ if the sum comprises all 
possible intermediate states and if in addition the isospin is conserved in the 
decay of resonances, one still has only two kinds of fragmentation functions: 



and 

^^-E^/--u.^/-i ■ (12) 



R 



d u 



These functions correspond uniquely to the standard favoured and unfavoured 
fragmentation functions and fulfil the relation (^-^) needed to obtain the 
identity (S). Thus we come to the conclusion that intermediate resonances 
do not violate Eq.(^) i.e. do not disturb the procedure of extraction of the 
d — u asymmetry. 

However, one cannot avoid completely the explicit treatment of interme- 
diate resonances when modelling pion spectra. The fragmentation functions 
(y) with the resonance components (|lO|) do not (!) obey the QCD evolu- 
tion equation; this is usually ignored in the current literature. It would be 
useful to separate out the direct fragmentation contribution to the fragmen- 
tation functions, which has a better chance of obeying the QCD evolution 
equations. However, to determine the fragmentation functions of the direct 
fragmentation one would need to perform a combined analysis of fragmen- 
tation into pions and into all other resonances having pionic decay channels. 
Such an involved experimental analysis has never yet been done. 



2.3 Choice of fragmentation functions 



In order to estimate the effect of nonpartonic components on the extraction 
of d — u asymmetry we need to fix fragmentation functions with which the 
main partonic term will be calculated. 

Modern parametrizations of fragmentation functions are fitted mostly 
to data from e~^e~ collisions. These analyses include leading or next-to- 



leading order QCD corrections (see for instance [15, 16, 17|). In contrast 
to e+e~ collisions the situation in ep scattering is much less developed: less 
experimental data, no QCD analysis. 

Let us see how the existing parametrizations of fragmentation functions 
behave in ep collisions. We start the quantitative estimations by comparing 
the existing fragmentation functions with the ^-distributions of charged pi- 
ous in ep scattering. In Fig.g we present ^ ( -^ — h '^^ ) data obtained long 
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Figure 2: Multiplicity distribution of charged pions -(-^j- 
Different sets of fragmentation functions are confronted with the Cornell JI^/ 
data with 2 < Q^ < 6 GeV^ and W ^ 3-4 GeV (left panel) and the EMC 
\19i] data with 20< <Q^> <71 GeV"^, depending on z (right panel). 



ago at Cornell jl8| ] (panel a) with kinematics similar to the HERMES exper- 
iment, and the data from EMC |l9| (panel b) with slightly higher Q^ . We 
show also the QPM predictions obtained with fragmentation functions from 
the fit to e"'"e~ data |]16|, from the fit to e^e~ and photoproduction data 
[P0| which include QCD corrections, and with fragmentation functions from 



the simple QPM fits to the ep data |21, 1£]. Surprisingly the "advanced" 



parametrizations give a much worse description of the data than simple ones 
do. However, simple parametrizations are limited to the relevant values of 

In principle the modern fragmentation functions used in e+e |16] were 
obtained including QCD corrections, i.e. beyond the naive quark-parton 
model. The correct formulae for the cross section in DIS calculated includ- 
ing QCD corrections are more complicated than QPM ones [^, ^. On the 
other hand the analysis of the HERMES experiment R] was performed based 
on simple QPM formulae, using Eq.(g). Therefore in order to compare with 
those results we have to stay at the QPM level too and ignore some inconsis- 
tency. Moreover the QCD evolution of fragmentation functions (included in 
the calculation) changes the pion multiplicity in agreement with the trend 
of experimental data and might create the main part of the Q^-dependence 
even when used with QPM formulae. The whole effect of the inconsistency 
should, however, be clarified in the future. 

In electron-positron scattering the number of negative and positive pi- 
ous produced is identical. This is not the case for ep scattering. Here quark 
distributions in the proton, isospin-asymmetric by their nature, allow us 
to distinguish between the favoured and unfavoured fragmentation which is 
very difficult, if not impossible, in e"^e~ scattering. As can be seen from 
Eq.(|8|) such a separation is essential for our analysis. In Fig.|3| we show z- 
distributions of negative (upper panel) and positive (lower panel) pions as 
measured by the EMC collaboration [19|. Different sets of fragmentation 
functions are confronted with the experimental data. A surprisingly poor 
description of the data is obtained with fragmentation functions from e+e~ 
scattering ||lq| . Not very good agreement of the Field- Feynman parametriza- 
tion is most probably due to a different Q^ here (^ 20 GeV^) than that 
where it was designed [21|. A correct QCD evolution should resolve this 
disagreement. As will be discussed below both QPM-parametrizations give 
reasonable representations of the ratio of unfavoured to favoured fragmen- 
tation functions, which is a more slowly QCD-evolving quantity. 

The ratio D^{z)/D^{z) directly enters formula (^ and, as our analysis 
shows, the measured d — u asymmetry is very sensitive to it. In Fig.Q we 
display the ratio obtained from the analysis of experimental data from EMC 
3[ and recently obtained by the HERMES collaboration at DESY 1 34 1 . The 



simple QPM parametrizations [£^, 19 1 provide a reasonable description of 



the data. In contrast the "advanced" fragmentation functions |1^, |20| fail 
again. Is it due to a different physics in e^e~ collisions than in DIS, or is it 
due to QCD corrections, or is it due to something else? In our opinion this 



10 



b N 



0.1 



0.01 
10 



b N 



0.1 




0.01 



— EMC parainetrization , v \ 

Field-Feynman \^^ \ 

- — Binnewies-Kniehl-Kramer '\ > 
Carlson-Wakely \\ 



0.0 



0.2 



0.4 



0.6 



0.8 1.0 

Z 



Figure 3: -(-^j— ) (upper panel) and -(-^j— ) (lower panel). The experi- 
mental data are taken from yjT^. The fragmentation functions used are the 
same as in the previous figure. 
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Figure 4: The ratio of unfavoured to favoured fragmentation functions as 
a function of z. The experimental data are from EMC |^/ and from a 
preliminary, unpublished HERMES analysis ^ 



is mainly due to the fact that the e~^e~ data are not sufficient to separate 
unambiguously the favoured and unfavoured fragmentation functions. A 
QCD analysis of DIS fragmentation functions is called for. 



The analysis above advocates the Field- Feynman parametrization [21| as 
the only good representation of the available ep data in the HERMES kine- 
matical region. This parametrization will be used in the following analysis. 



3 Nonpartonic components 

For small Q^, as for the HERMES experiment, some mechanisms of non- 
partonic origin (see e.g. Fig.|5|) may become important. For instance, the 
virtual photon can interact with the nucleon via its intermediate hadronic 
state. Such a mechanism is usually described within the vector dominance 
model (VDM). The photon could also fluctuate into a pair of pions, where 
both or one of them interact with the nucleon. In addition some exclusive 
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Figure 5: Nonpartonic mechanisms of pion production taken into account 
in this work: (a) VDM contribution, (b) spectator mechanism, (c) elastic 
production of the p" meson and its decay. 

processes can produce pions directly or as decay products of heavier mesons. 
To our best knowledge none of such processes has been investigated in 
the literature. Their influence on the extracted d—u asymmetry also remains 
unknown. We shall discuss processes shown in Fig.g one by one. 



3.1 Central VDM contribution 

Let us start from the VDM component (see Fig.|5|a). It was shown that 
in the inclusive DIS incorporation of the VDM contribution and related 
modification of the partonic component help to understand the behaviour 
of structure functions F^ and F2 at small Q^ pT| , p2| . This model was 
confirmed by a recent analysis of the Q^ dependence of the world data for 
the structure function difference F2 — F2 ||5[ . 
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The model for inclusive structure functions [31, 32] can be generalized 
to semi-inclusive production of pions: 

^ v'yv {Q +Mv) 

The second sum above runs over vector mesons V = p^,uj,4> and Oy is 
a correction factor which takes into account finite fluctuation times of the 
virtual photon into vector mesons for large x |32|. 



The inclusive cross section for pion production in vector meson {p^,u!, <p) 
scattering off the proton and neutron is not known experimentally. There is 
no model in the literature that one can trust quantitatively but in analogy to 
the total p^N cross section the cross section p^N — > ir X can be estimated 
as: 

a{p°p -^ TT^X) ^ 1/2 [a{TT+p -^ TT^X) + a{7r-p -^ vr^X)] 

a{p^n -^ TT^X) « 1/2 [cr(7r+n -^ tt^X) + a{TT~n -^ n^X)] . (14) 



a{p°n - 


■^vr+X) = 


-- <r{p'p- 


-^TT~X) 


a{p^n - 


■^7r~X) = 


= C7{p%- 


■^TT+X) 



Experimental data from the ABBCCHW collaboration [|36[ at pj^jj = 8, 
16 GeV correspond approximately to the range of the HERMES experiment 
[^. Unfortunately as it often happens in high-energy physics there is only 
data for proton targets. Using isospin symmetry for hadronic reactions one 
can obtain corresponding cross sections on the neutron from those on the 
proton by assuming 

(15) 
(16) 

These relations hold not only for the total cross sections but also for differ- 
ential ones independently of energy. From the most complete data at pf^^^ = 
16 GeV [H] (W = 5.56 GeV) we get 

1/2 [a{7T+p -^ TT+X) + a{7r-p -^ tt+X) ] = 38.65 ± 0.29 mb , 
1/2 [a(7r+p -^ TT-X) + a{-K-p -^ vr^X) ] = 31.80 ± 0.22 mb ; 
clearly different values. Although the bulk of the difference comes from the 
target fragmentation region, which we are not interested in, in the beam 
fragmentation region it was also found (t[tt'^p -^ tt~X) = 14.8 mb and 
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a{'K~p -^ TT^X) = 19.0 mb [^] with almost equal cross sections for beam- 
like pions. 

The situation in the semi-inclusive case is more compHcated than for 
total cross sections. The experimental spectra for vr^p -^ n^X contain 
components due to peripheral processes, which are, in general, specific, dif- 
ferent for different reactions. Peripheral processes from the vr+p -^ tt~^X and 
-K^p -^ -K^X reactions do not contribute to the p^p -^ tt X reaction and 
should be eliminated; only nondiffractive components of the irp -^ irX re- 
actions should be taken into account when modelling p'^-induced reactions. [| 
This requires physically motivated parametrization of the vrA^ -^ nX data. 

Following these arguments we have parametrized the experimental dif- 
ferential cross sections for four different reactions it p -^ tt X from |36| as 
a sum of central and peripheral components 

da da'"''' daP^"- 

+ :^-r^' (17) 



dxpdpj^ dxpdpj^ dxpdp^ 

where xp is the well known Feynman variable. Details of this analysis will 
be presented elsewhere, a short sketch is given in Appendix ^. Because the 
CM-energy of the ABBCCHW collaboration is very similar to that of the 
HERMES experiment, we believe that in the range of energy relevant for 
this experiment the functional form given in the appendix is suitable. 



Finally, our analysis of experimental data [36| combined with the as- 
sumption of isospin symmetry (^6|) shows that for the nondiffractive com- 
ponents still 



a{p"p -^ TT^X) / a{p"n -^ ir^X) . (18) 



In Fig.^ we show 



„o 



^^^ _ ajp^p ^ TT^X) - ajp^n ^ vr^X) 
V ^^pOp ^ ^±x) + a{pOn -^ 7r±X) ^ ' 

for both positive and negative pions. An identical value of the asymmetry is 
expected for the co beam and a similar one for the 4> beam. This automati- 
cally means practically the same result for photon (real or virtual) induced 
reactions which proceed via hadronic intermediate states. We obtain rather 
large asymmetries, larger than for total photoproduction cross sections on 
the proton and neutron. 

^ Some peripheral processes specific for the p" beam will be included explicitly in 
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Figure 6: (5pn~"^ calculated from Eq.^^ as a function of Feynman xp. The 
pion+proton cross sections were taken from the parametrization described in 
Appendix 0. 



Different cross sections on proton and neutron mean that the VDM con- 
tribution modifies the r.h.s. of Eq.(|8|). In Fig.|^ we show a modification of the 
measured quantity ^5f due to the central VDM component. In the HER- 
MES experiment [^ both Q^ and W vary with Bjorken-x, but the change of 
energy is considerably smaller. In the present calculation the photon-proton 
CM energy was fixed at the average HERMES value W = 5.0 GeV^. In the 
panel (a) we show a result of a calculation where the central VDM com- 
ponent discussed in this section is simply added to the main fragmentation 
component. In the panel (b) the fragmentation component in addition was 
rescaled by the factor ^2,^2 (see Eq.(0)) which is required for consistency 
of inclusive and semi-inclusive structure functions. The solid line in both 
panels represents ;j^ obtained directly from the parton distributions |64]. 
As can be seen from the figure the r.h.s. of Eq.(|8|) clearly deviates from the 
assumed partonic outcome. The effect is surprisingly large, especially for 
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Figure 7: The true (solid) and the modified by the central VDM contribution 
:[^, calculated according to the l.h.s. and r.h.s. of Eq.^, respectively, as a 
function of Bjorken-x for different values of z and typical HERMES W = 5 
GeV. The result in panel (a) is obtained by a simple addition of the fragmen- 
tation and the central VDM components (Qq = in Eq.^T^)). The result 
in panel (b) is obtained with the rescaled (Qq =0.8 GeV"^) fragmentation 
component, as described in the text. 
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small x^ Thus, the quark flavour asymmetry extracted from semi-inclusive 
experiments in the simple QPM approach seems to be highly overestimated 
if the VDM contribution is neglected. 

The VDM effect discussed in this section is not completely new. A 
similar effect of the hadronic structure of the photon on the difference of 
semi-inclusive cross sections o\yp_^^+x — cr^p^^-x was already noticed long 
ago in real photoproduction |^2|. Although in DIS the effect is smaller, it, 
however, strongly modifies the measured d — u asymmetry. 

3.2 Spectator mechanism 

In both partonic and central VDM mechanisms the virtual photon is to- 
tally absorbed and pions are produced in a complex process involving many 
degrees of freedom. Such pions are then preferentially emitted at not very 
large values of z. The peripheral processes are not included either in the 
partonic or the central VDM component considered above and are expected 
to be important in the region of large z where these processes disappear. 
Let us consider first the spectator mechanism depicted in Fig.^D. To our 
knowledge such a mechanism has been never discussed in the literature for 
photon induced reactions. We begin with the case of real photoproduction 



where we can apply a technique from [37|. For virtual photons with g^ < 
the formalism is not well developed. 

The cross section for the spectator pions vr^ in real photoproduction can 
be expressed as a product of the distribution of pions in the photon {f^±/^) 
and the off-shell total cross section for scattering of tt off the proton or 
neutron: 

-^^L./,{z).a^:^{il-z)s). (20) 

At small vrA^ energies St^n ~ {1—z)s^*n relevant for the spectator mechanism 
there can be a difference between cr^^^^ and cr^^^^. Together with isospin 
symmetry of hadronic reactions this would lead to different N^ and N^ 
i.e. would modify the identity (^). To represent the total cross sections of vr^ 
scattering off nucleons we use cubic interpolation of the world experimental 
data [^1. 



'^ We wish to remind the reader here that in the HERMES experiment the photon 
virtuahty for small x is only of the order of 1 - 2 GeV^ . 
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Figure 8: The invariant structure function /(xp) for the 'yp —^ n X reaction 



and E, 



7 



9.3 GeV. 



The distribution function /.^±/^ can be calculated using a technique sim- 



ilar to that outlined in |37]: 
f^±l^{z,ki) 



9em 



647r2 zil 



1 



2kl 



[0-M2Jz,fei)]2 



\F{zM 



(21) 



where gem is the electromagnetic 7 — > vr+vr" coupling constant, Mt^tt is 
the invariant mass of the two-pion system, F{z, k"^) is a vertex form factor 
which accounts for the finite size of particles involved and off-shell effects. 



For other technical details see [l37 |. 

In Fig.^ we show the result of our calculation (dashed line) with cut- 
off mass 1.5 GeV as in Ref.|37] for s^" = 4.28 GeV corresponding to the 
experimental data from SLAC [^] with E^ = 9.3 GeV. The results are pre- 
sented in terms of the invariant single particle structure function [^^ (see 
also Eq.(^3|) in Appendix^. The energy in the pion-proton subsystem de- 
creases with increasing xp and we observe fluctuations due to s-channel ttN 
resonances. The peaks would be even more pronounced at smaller photon 
energies and would disappear completely at larger photon energies. 
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The contribution of the spectator mechanism calculated with reasonable 
cut-off masses in the vertex form factor is smaller than experimental data. 
A dominant fraction of pions produced at large xp appears to be given 
by the mechanism of elastic p^-meson production and its subsequent decay 
(discussed in detail in section 3.3). By the short-dashed line we present this 



contribution corresponding to the cross section calculated as: 

^'"'^^^ = -(7P -A) •/£-:;' (^), (22) 

where the explicit functional form of the decay function f^^^ (-z) can be 
found in Appendix p. In the calculation we have used the experimental 
value (t(7p -^ p^p) = 13 /ib, relevant for E^ = 9.3 GeV. The sum of both 
contributions (solid line) seems to be consistent with the large-Xi;' part of 
the pion spectrum. Thus in real photoproduction the spectator mechanism 
becomes non-negligible only at large xp (or z). 

In virtual photoproduction the situation is more complicated. The abso- 
lute normalization of /tt* /-y* should depend on the virtuality of the photon. 
However, it cannot be calculated from first principles as it involves a ver- 
tex form factor where more than one particle is off mass shell. Empirically 
such cases are strongly damped pO| . Also a naive use of a typical light-cone 
parametrization of the vertex form factors |^3l with the mass of the parent 
particle replaced by the virtuality of the photon leads to a strong suppres- 
sion in comparison to real photoproduction. How strong this suppression is 
in comparison to the suppression of the elastic-/?'^ production mechanism is 
not clear. For photon virtualities of the order of Q^ ~ 2-4 GeV^ we find that 
the elastic-p*^ contribution, discussed in detail in the next section, together 
with the partonic component, totally accounts for the cross sections at large 
z, leaving practically no room for the spectator mechanism (see Fig.|lq). 

In summary , the spectator mechanism, while potentially important in 
real photoproduction is most probably negligible in DIS. 

3.3 Exclusive p meson production 

The exclusive meson production 'y*N ^ MN' is one more mechanism not 
included in the fragmentation formalism (Eq.(Q)) and may also modify the 
extraction of the d — u asymmetry. The pion exclusive channels (M = tt)q 

^ These are not included in the spectator mechanism discussed above where the final 
state X — N was not taken into account. 
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contribute at z ^ 1, i.e. outside of the range of the HERMES kinematics 
and will be ignored in the following discussion. In contrast, the pions from 
the decays of light vector mesons may be important in the context of the 
d — u asymmetry from semi-inclusive pion production. The production of 
p mesons (M = p) seems to be of particular importance. Firstly, the p^N 
channel is known to be the dominant exclusive channel in 7*A^ scattering. 
Secondly, because p^ decays predominantly into two pions this will produce 
pions with <z>~ 2- ^ detailed calculation (see Appendix p[) shows that 
the dispersion of the decay pion z-distribution is large and therefore this 
effect has a chance of being observed at large z where the hadronization 
rate is already much smaller. The next potentially important mechanism is 
the uj meson (M = a;) production. However, the dominant uj meson decay 
channel is the three-body system ■k^-k~-k^ i.e. it is expected to contribute 
to the inclusive pion distribution at < z > ~ g, i.e. in the region where the 
hadronization rate is large. Moreover the cross section for the uj channel 
is smaller than for the p" channel. Below we shall consider the p^ channel 
only, which is probably the most important. 

The elastic p^'-production contribution (diagram (c) in Figj5|) to the semi- 
inclusive structure function (^Jl3|) can be written formally as 

J-f '^'(x, Q\ ^) = ^ • ^7*^-P°^(^' Q^) ■ fp^^^A^) ' (23) 

The decay function /pO_»^ can be easily calculated (see Appendix ^. For 
not too high energies, as for the HERMES experiment one may expect 
(t(7*p -^ p^p) 7^ a{'y*n — > p^n) which would modify the d — u asymme- 
try. At high energy the pomeron-exchange (two-gluon exchange) mecha- 
nism dominates and one may expect cr['y*p -^ p^p) = a{'~f*n -^ p^n). At 
low energy the exchange of subleading reggeons (quark exchange) could lead 
to a{'y*p -^ p^p) / a{'y*n — > p^n) due to isovector contributions. In real 
photoproduction the isovector amplitude is known to be rather small J43U . 
In DIS the situation may be quite different. Assuming that the production 
is hard, i.e. of perturbative nature, the longitudinal p^ is predicted to be 



dominated by the quark exchange mechanism at low photon energies |44]. 
The HERMES 7*p energy corresponds precisely to the maximum of the 
p^ production in the hard quark-exchange exclusive reaction 7*p -^ p^p 



[44, 45, |4g, 47 1 . Although there are no experimental data in this region, the 



data from EMC |61[1, NMC [|5| and E665 |62|] collaborations in the close 



neighbourhood seem to be in rough agreement with the calculation for ai 



[44, 45 1 . Different quark distributions in the proton and neutron lead in this 
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approach to different p^ production cross sections for proton and neutron 
targets, which obviously leads to a different production rate of charged pions 
in reactions on the proton and neutron. This in turn modifies the right hand 
side of Eq. ^ and the subsequent conclusions on the d — u asymmetry. 
One could also try to understand elastic meson production within the 



Regge phenomenology ||49|, 50 1. It is not obvious a priori what is the kinemat- 
ical range of applicability of either the quark exchange approach or Regge 
phenomenology. Below we shall investigate elastic p^ production on the pro- 
ton and neutron using both these approaches: 

a) the Regge approach, and 

b) a QCD inspired quark-exchange model. 

3.3.1 Regge approach 

The cross section for neutral p meson electroproduction in the HERMES 
kinematics is not known experimentally. While for the proton target there 
are data in slightly different kinematical regions of x and Q^ |5^, there is 
almost no data for the neutron target. Only in one work ||5J] was the p^ 
production studied simultaneously for the proton and deuteron data. The 
X- and Q^-integrated result obtained there 



Tincoheritd ^ pV) = ((0-7 - 0.8) ± 0.2) • a(7> ^ A) 

does not exclude the difference between the proton and neutron target. If 
the nuclear effects are completely ignored this leads to 

o-(7*n -^ p°n) = ((0.4 - 0.6) ± 0.4) • ai'y*p -^ p^p). 

However, in the present analysis we are interested in x- and Q^-dependent 
cross sections. The only electroproduction data on the deuteron target with 
well defined kinematics were published in [£^]. Such a limited set of the 
deuteron data does not allow us to determine the cross section on the neutron 
target for x and Q^ in the kinematical region we need. A possible way 
out would be to parametrize the proton data with a suitable, physically 
motivated parametrization in a possibly broad kinematical range (there are 
rich data around the kinematical region of the HERMES experiment, see 
Fig.|9|) and use theoretical arguments and/or experimental data for other 
reactions to determine the neutron cross sections. 
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In the present paper we shall parametrize the existing experimental data 
for exclusive p^ production by means of the following simple Regge-inspired 
reaction amplitude, similar to that in Ref.p9|, 



^Xz'>::{i*N^p'N-t) 



i ■ Cjp{t) ( - 



± 




Cis{t) - 



with " + " and " — " in front of the isovector {IV) contribution for pro- 
ton and neutron, respectively. The pomeron contribution is marked by IP 
and isoscalar reggeon contribution by IS. The following normalization is 
assumed: 



da 
H 



1 



{2sN + l)Nx 



■E 

A-,* Xy 



A 



^AT 



7 

-Ajv 



it) 



(25) 



-1/2 



(24) 



where sn is the spin of the nucleon and N\ is the number of active helicity 
states of the virtual photon. 

In the following we are interested in relatively low 7*A^-energies where in 
principle the pion-exchange mechanism could be important too. At low en- 
ergies, just above resonances, the pion-exchange mechanism is known to be 
the dominant mechanism for photoproduction of u; mesons ||5l|, 52 1. It can be 
shown that due to the helicity structure of its amplitude the pion-exchange 
contribution does not interfere with the diffractive contribution as far as the 
spin-averaged cross section is considered, that is the pion-exchange contri- 
bution can be added incoherently in the cross section. Because r„o^^o , ^ ^ 



^7r0+7 the relevant coupling constant f^o^ 



967r 



— ^ T 



Fp^^ is rather 



small. In order to estimate the corresponding cross section one has to make 
some plausible estimation for the vertex form factors. Assuming the form 
factors which lead to a good description of the oj-photoproduction data, 
whilst neglecting other mechanisms, provides a reasonable upper estimate 
of the pion-exchange contribution for p^ production. The pion-exchange 
contribution estimated in this way can most probably be neglected in the 
kinematical region considered here. It is interesting to notice that unlike 
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for the neutral p meson production the pion-exchange contribution cannot 
be neglected for charged p meson production due to lack of the dominant 
isoscalar contributions. The discussion above further justifies the simple 



Ansatz used in (24). 

In practical application we assume the same t-dependence of C]p, Cjs 
and Cjv and take A = nip. The total 'y*N ^ p^N cross section can be 
obtained as the integral 

a{W,Q^)= ^{t)dt. (26) 

There are many approximate or even incorrect formulas for the upper and 
lower integration limits in the literature. It is, however, essential to use 
correct formulas for small W . 

The free parameters in Eq.(24) i.e. e^p, Cjp and Cjs + Cjv have been 



fitted to the existing experimental data for the angle-integrated cross section 
for p^ production on hydrogen [^]. The slope parameter B in the exponen- 
tial i-distribution was fixed at -B = 6 GeV~^ which is known experimentally.^] 
In order to avoid poorly understood contributions of baryonic resonances we 
have limited our fit to VF > 3 GeV. The simple form of the amplitude above 
is obviously not adequate for large Q^ where genuine hard QCD processes 
take place. Consequently we have included in our fit only experimental 
data points with Q^ < 10 GeV^. In Fig.^, together with all available p^ 
electroproduction data, we have shown these kinematical boundaries. The 
large filled circles denote the kinematical loci where the HERMES analysis 
of charged-pion semi-inclusive data was made. With the above cuts we get 
from the fit: Cjp = 1.57 pb^/^GeV'^, Cis + Civ = 6.33 pb^/^GeV'^ and 
ep = 0.102. The quality of the fit is shown in Fig.lC. The corresponding x^ 
= 2.34. 

To separate the sum Cjg + C/y into the isoscalar and isovector parts one 
needs more information. The size of the isovector a2-exchange contribution 
was estimated long ago for total photoproduction cross section (see for in- 
stance |4^). It was found empirically that the total photoabsorption cross 
section on the proton and neutron can be parametrized at low energy as 

otot = Cjp + {Cf ± Ca, )E-]liB ■ (27) 



^It is impossible to separate the isoscalar and isovector reggeon contributions from the 
fit to the proton data only due to their identical energy dependence. 

® The results are rather stable against a small variation of B in the range B = 6 ± 2 
GeV-2. 
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Figure 9: Experimental data points for exclusive p^ production on the proton 
and deuteron (open squares). For illustration the HERMES kinematics is 
shown in addition by large solid circles. The thick solid lines show the limits 
for the Reg ge- inspired fit. 

In our parametrization of the j*p -^ p^p data the energy dependence of 
the pomeron- and reggeon-exchange contributions is slightly different (see 
Eq.(p^)). Extending the validity of Regge phenomenology to both real and 
virtual photons we can write somewhat schematically the amplitude 



A{YN -^ 7*iV) 
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02 



(28) 



for Compton scattering, 
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Figure 10: The cross section for j*p — > p^p as a function of center- of-mass 
energy for selected values of photon virtuality. The solid line is obtained 
from the VDM-Regge inspired fit. The dased line shows the cross section on 



the neutron target. The solid squares represent the NMC nuclear data [551. 
Please note that excepting Q^ = 0.47, all other curves and experimental 
points are resettled by 5, 5^ , 5^ and 5^. 
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^(7*iV -^ p^N) 



m7 



mf^ + Q^ 



for exclusive p^ photoproduction and 



A{j*N -^ ujN) 



mt 



m. 



+ g2 



—]P+—f±—a2 
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(29) 



(30) 



for exclusive oj photoproduction. We have used the values of 7p0 and 7^^ 



from 



and put nip 



rriu 



my. Please note that iP, / and 02 cor- 



responding to the reggeon-exchange amplitudes on the hadronic level are 
the same in all these reactions. Different factors in front of these hadronic 
amplitudes give different strength of each contribution in different reac- 
tions. We have adjusted the relative strength of the a2-contribution com- 
pared to the /-contribution in the Compton scattering amplitude ( p8| ) to 
reproduce the empirical low-energy parametrization (p7|) for atotiip) and 
atotim). In Fig.|T] we compare 5^^^'°'' = ZiZ^i^.TZf:!^, given by the ob- 



tained parametrization with that from the empirical fit. Shown is the band 
due to the uncertainties of parameters from P^ j and the best representation 
of the empirical formula by our Regge parametrization (p4[). Although we 
reproduce the empirical fit rather well, in our case there is a different rela- 
tive strength of the / and 02 contributions to the Compton amplitude. This 
difference is caused by the energy dependence of the pomeron exchange as 
opposed to the constant assumed in the empirical fit P^ ]. 

Having fixed parameters in (p8|) we can calculate the corresponding 



proton- neutron asymmetry for p^ and uj production: 5pn, 5p„, which are 



also shown in Fig. 11. We have shown in addition experimental results for p^ 



[p6| and w [57| photoproduction which are consistent with our parametriza- 
tion. While the asymmetry of the cross sections for the p^ production is 
similar to the Compton case, the asymmetry for the lo production is con- 
siderably larger. [| The latter may also be important in the context of J — n 
asymmetry. However, the absolute normalization of the corresponding cross 
section for the j*N ^ ujN reaction is not well known, at least in the region 
of the HERMES kinematics. It is expected to be considerably smaller than 
for the p^ production. 



^ We have neglected here the pion-exchange mechanism which would decrease 5^ 
energies W^ < 5 GeV. 
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Figure 11: The asymmetry Spn as a function of center- of-mass energy. The 
thick solid line is obtained from the Regge parametrization adjusted to re- 



produce the empirical fit for total photoproduction (27) represented by the 
hatched band. The dashed and dotted lines correspond to a similar asymme- 
try for p^ and lu production, repectively, and were obtained from our Regge 
parametrization. The experimental points are taken from pqj and Wi 



Although the difference of the cross sections for exclusive p^ production 
on the proton and neutron targets is small, the effect of this mechanism on 
the d — u asymmetry is not negligible at all. In Fig.^, we show the corre- 
sponding modification of the quantity ^^^ in the same way as before for the 
central VDM contribution. This modification may even be underestimated, 
as it is based on the Regge-inspired parametrization of the cross section for 
the 7*n — > p^n reaction which overestimates the p^ production on nuclei. 
This can be seen in FigJin where the NMC experimental points lie below 
the Regge-parametrization. 
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Figure 12: T/ie trite (solid) and the modified by the exclusive p^ production 
{d—u)/{u—d) as a function of Bjorken-x for different z and typical HERMES 
W = 5 GeV. As in the central VDM case, in panel (a) Qo = ^ (standard 
partonic component) and in panel (b) Qq = 0.8 GeV'^ (rescaled partonic 
component). 

3.3.2 QCD approach 

For sufficiently large Q^ ^ the elastic vector meson production can be cal- 
culated perturbatively in the formalism of off-forward parton distributions 
(OFPD's) [Q, ^^. It was argued that only the cross section for longitu- 
dinally polarized photons, where end-point contributions are suppressed, 
can be calculated reliably |5^. The cross section for the exclusive reaction 
'y'^ + N ^ Vl + N can be calculated in the standard way as 



da 



■y+N^V+N 
LL 

dt 



1 1 



16tts^2 ^ 






(31) 



It is not clear at present how large the virtuality should actually be for the applicability 
of this formalism. 
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where Xn and A'^ are helicities of the incoming and outgoing nucleons, 
respectively. The two zeros in the upper index row of the matrix element 
correspond to longitudinal photons and helicity of the produced vector 
meson. The amplitude of the two-quark exchange mechanism for vector 
meson production was calculated for the first time in p4, ^ ]. The total 
longitudinal cross section can be obtained by integrating (3l|) over t in the 
kinematically allowed interval. 



In the formalism first proposed by Ji |60|, neglecting transverse momenta 
of quarks in the nucleons and in the vector meson, the leading order ampli- 
tude reads |45] 






9QJo z(1 - z) 2 



dx 



1 



+ 



^ + ie X + S, 



le 



(4vra,) H^{x, ^, t) N{p', X'^)j ■ nN{p, Xn) , (32) 



where ^v{z) is the distribution amplitude and H^{x,^,t) is a generalized 
function related to so-called skewed quark distributions in the nucleon. For 
the electroproduction of p^ mesons we are interested in here one gets 



HP^°{x,^,t) = ^ 



1 

72 



V/^(x,e,t) + V/^(x,e,t) 



(33) 



The functions iJ"/^(x,C,t) = UN{x,^,t) and H'^/^{x,^,t) = dN{x,^,t) are 
the non-diagonal, off- forward quark distributions. In this subsection we shall 
concentrate on the relative magnitude of the cross sections for p^ production 
off the neutron and proton. Therefore the approximation relying on the 
replacement of .^ ^ x (i.e. using familiar diagonal quark distributions) 
seems sufficient for our present purpose. Consequently we shall take 



for X > and 



H^/^ix,^,t) 
^'^/^(x,e,t) 



H''/^{x,i,t) 



un{x) ■ D{t) 
dN{x)-D{t) 



-UNix) ■ D{t) 
-djv(x) • D{t) 



(34) 



(35) 



for X < 0. The above Ansatz assumes factorization of x and t dependences. 
Thus the whole t-dependence will be contained in one universal function 
D{t), common for all flavours. We shall try exponential and dipole form 
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factors which provide a good representation of experimental data for exclu- 
sive /?" production. The factorized form has the advantage that the total 
longitudinal cross section can be obtained analytically. 

In the present analysis we have neglected the tensor magnetic-type E- 
terms (see [^]) ^ which may be expected to be important only at large t 
and lead therefore to a rather small contribution to the total cross section. 

The integral in the amplitude given in ( |32|) can be calculated in the 
standard way by splitting the integral into a real principle value and an 
imaginary (5-function. 

The cross section asymmetry defined as 






pn -, * p — ^ pO 10 --y * r) — * rS' n ^ ' 



calculated according to (|3l|) and (|3^) is shown in Fig.|^ as a function of 
photon-nucleon center-of-mass energy W for Q^ = 4 GeV'^, typical for the 
HERMES kinematics. As in previous calculations the quark distributions in 
Eqs.p^) and (35) were taken from p3. If only valence quark distributions 



are taken into account there is a relatively large asymmetry between the 
scattering off the proton and neutron targets. The inclusion of sea quarks 
decreases the asymmetry, which vanishes completely at large energy (small 
X for a fixed Q"^). At the energy of the HERMES experiment PF ~ 5 GeV 
there is about 10 percent asymmetry i.e. about a factor 2 more than in the 
Regge approach. 



In Fig. 14 we compare the asymmetry obtained within the Regge ap- 
proach discussed in the previous section and within the perturbative for- 
malism discussed here for fixed W (left panel) and fixed x (right panel) as 
a function of Q^. There is a substantial difference between the Regge ap- 
proach, where W is the variable relevant for the proton and neutron asym- 
metry, and the QCD approach, where it is rather Bjorken-x which controls 
the asymmetry. In the Regge approach, if W is fixed the asymmetry is 
practically independent of Q^ and varies strongly for fixed Bjorken-x. In 
contrast, in the perturbative approach for fixed Bjorken-x the asymmetry 



only weakly depends on the photon virtuality, as can be seen from Fig. 14b. 
If the center-of-mass energy is fixed instead, the Q^ dependence of the ra- 



tio is much stronger. As can be seen from Fig. 14 the asymmetry between 



proton and neutron target becomes larger for larger photon virtuality. 



As far as we know these terms were never estimated in the literature. 
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Figure 13: (5p„ for quark- exchange mechanism as a function of center-of- 
mass energy. The solid line includes both valence and sea quarks. For com- 
pleteness we show also the asymmetry for valence (dashed line) and sea 
(dotted) quarks exclusively. 



We have calculated only longitudinal cross sections which dominate at 
large Q^. At small Q^ the transverse cross section becomes equally im- 
portant. Although it is not possible to make a rigorous calculation for the 
transverse cross sections, it is natural to expect 6T^ and 6p^ to be simi- 
lar. Therefore we expect that the asymmetry for longitudinal cross sections 
should be a reasonable estimate of the asymmetry of the total (longitudi- 
nal-|-transverse) cross sections. 

The description of the experimental data for exclusive p^ meson produc- 
tion by means of the hard mechanism (not discussed here) is not as good as 
by means of the Regge approach. The absolute normalization of the cross 
section depends on transverse momentum distributions of quarks in the nu- 



cleon and in the produced p meson [45] which are not fully understood at 
present. Therefore we shall not calculate here the corresponding modifica- 
tion of the measured {d — u)/{u — d). Such an analysis requires first of all 
a good description of the absolute value of the cross section. We expect, 
however, at least as big a modification as in the Regge case. 
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Figure 14: 5p„ as a junction of Q^ in the Regge (dashed) and QCD inspired 
(solid) approaches for: different fixed W = 3, 5, 7 GeV (left panel) and 
different fixed Bjorken-x = 0.05, 0.1, 0.2, O.4 (right panel). 



3.3.3 Charged p mesons 

Above we have considered only neutral p mesons. Charged p mesons can 
also be a source of charged pions due to their decay mode p^ —>■ tt^tt^. 
Experimentally the cross section for exclusive charged p mesons is much 
less well known than that for the neutral p mesons. It can be estimated 
within the QCD formalism of the OFPDs approach as that sketched for 



neutral p mesons using symmetry relations for the matrix elements |47]. 
Because the cross section for exclusively produced charged mesons depends 
on quark distributions in the nucleon differently than in the QPM formula 
(^) the contribution of charged p mesons will certainly modify the d — u 
extracted by means of (|^). Because the result depends on the magnitude 
of the charged-meson production which is rather difficult to predict in the 
QCD-type calculations (off-diagonal effects, the choice of the scale of the 
running coupling constant, inclusion of transverse momenta) we shall leave 
the problem for a separate more refined analysis. 
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4 Comments on nuclear effects in the deuteron 

So far we have followed the HERMES collaboration and neglected all nuclear 
effects in the deuteron i.e. assumed that 

a{-f*d -^ vr^) = a{-f*p -^ vr^) + a{j*n -^ vr^) . (37) 

The theory of nuclear effects in semi-inclusive processes is less developed 
than in the inclusive case. Let us consider a simple example of an x- 
independent relative nuclear effect of size k, universal for produced vr"*" and 
vr", 

a{-f*d -^ 7r±) = (1 - k) {(t{j*p -^ vr^) + a{-f*n -^ tt^)) . (38) 

Then the semi-inclusive cross sections for pion production on the neutron 
extracted from the deuteron target data are 

= (1 - K)o-(7*n -^ -K^) - Ka{j*p -^ TT^) (39) 

= cr(7*n — > vr ) — k (cr(7*p — > vr ) -|- a{'y*n ^ vr )) 

i.e. biased by the assumed nuclear effect k in the deuteron. 

Thus the differences o"(7*p -^ vr ) — (T(7*n -^ vr ) needed in Eq.(|8|) are 
replaced by 

= {l + K)a{j*p -^ 7r±) - (1 - K)a{-f*n -^ 7r±) (40) 

= (j{'j*p —>■ vr^) — a{'y*n — > vr^) -|- k {a{'y*p —>■ vr^) + a{'y*n -^ vr^)) 



In Fig.|T5|we show the nuclear effects on ^ for k = 0.02, 0.01, 0.0, -0.01, 
-0.02, i.e. in the range known from inclusive DIS. These effects are indepen- 
dent of z by assumption (pH). Following the inclusive case, for small values 
of Bjorken j; < 0.1 a shadowing, i.e. k > is expected which means that the 
asymmetry obtained when neglecting nuclear effects is underestimated (see 



Fig.|l5|) 



The shadowing leads to an effect opposite to that for the resolved photon 
component discussed earlier in this paper. For somewhat larger x an anti- 
shadowing due to excess pions is not excluded. For still larger x a nuclear 
binding and Fermi motion corrections come into play. 

Summarizing, we have shown that even small nuclear effects, of the or- 
der of just a few percent, lead to considerable consequences for the d — u 
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Figure 15: The true (solid) and the modified by the nuclear effects (d — 
u)/{u — d) as a function of Bjorken-x for W = 5 GeV and different values 

of K. 

asymmetry. Nuclear effects are expected to be x and Q^ dependent. In the 
present analysis we have shown only a band of uncertainties due to nuclear 
effects. A more precise determination of the x or Q^ dependence requires 
a more microscopic calculation which goes beyond the scope of the present 
paper. This is, however, necessary if the d — u asymmetry is to be extracted 
from semi-inclusive data. 

5 Consequences for the HERMES experiment 

Having discussed each of the nonpartonic effects separately we will now 
attempt to combine them and try to understand their net effect on the 
measured d — u asymmetry, and the consequences of this for the HERMES 
data 0. 

Before a numerical estimation of the total nonpartonic effect we would 
like to discuss briefly a subtle problem. At first sight it seems that fragmen- 



33 



biN 



bl N 
xJro 



piarc part- approach 
VDM 

elastic p 
part.+VDM + el.p° 
part." +VDM + el,p 



10 



Ir- 



0.1 ^ 




OOOOO Cornell data 
(Q') = 2.8GeV'~ 




0.01 r 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 16: Multiplicity distribution of the charged pions. Contributions of 
different mechanisms are shown separately. The total effect calculated in two 
different ways discussed in the text is also presented. 



tation functions fitted to experimental data effectively contain nonpartonic 
effects. This is most probably not true as nonpartonic components are 
higher-twist effects i.e. are strongly Q^-dependent in contrast to the leading 
twist. 

In section p.3| we have selected the fragmentation functions which de- 
scribe the pionic yields well. There is a danger a priori that explicit inclu- 
sion of the nonpartonic effects discussed in the present paper may worsen 
the description of pionic spectra. In Fig.lC we show the multiplicity of 



charged pions as a function of z for <Q^ > = 2.8 GeV^ (compare with the 



left panel of Figj2D. Together with the corresponding Cornell data [18| we 
show the contributions due to different mechanisms separately. The contri- 
butions of nonpartonic mechanisms are considerably smaller than the main 
partonic contribution, however not negligible. If we add all of them together 
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we obtain the multiplicity (thick short-dashed line) over the experimental 
points. Another way to incorporate the partonic and nonpartonic compo- 
nents was proposed in Ref.|B^] for the inclusive case. As can be seen from 



Fig. 16, extension of this model to the semi-inclusive case, i.e. rescaling of the 
partonic component as in Eg. (]l3|) , (thick long-dashed line) provides a very 
good description of the experimental multiplicities. This approach treats 
all contributions explicitly, which seems more consistent than the approach 
mentioned above including them all effectively into fragmentation functions. 



It is also consistent with the inclusive structure function model |32|. 

In Fig. 17 we show ^^5^ as a function of z for different bins of Bjorken- 



X (the averaged value of x is given in the figure). The data points are 
taken from 0. The filled squares correspond to the data averaged over z. 
The QPM (z-independent) prediction with leading order quark distributions 
from |6J] is shown for reference by the solid line. The results of the calcula- 
tion including partonic, VDM and elastic (P contributions are given by the 
dashed line (with the QPM contribution calculated in the standard way) 
and long-dashed line (with the QPM contribution modified as in Eq.(|l3|)). 
As in the previous calculations the Field-Feynman fragmentation functions 
PH were used here. 

As can be seen from the figure there is a significant deviation of the 
"measured" asymmetry from the "real" one (the difference between the 
dashed and solid lines), especially for small values of Bjorken-x. The effect 
is bigger for the "modified-QPM" approach i.e. for the more consistent one. 
The nuclear effects most probably will introduce further deviation which is, 
however, difficult to estimate numerically. A significant z-dependence casts 
doubts on the averaging in z at least in the whole range from to 1. How- 
ever, as seen from the figure the experimental statistics and the z-range of 
the HERMES experiment do not allow this dependence to be identified. 

It is worth noting that the modification of the QPM result depends 
strongly on the fragmentation function used to calculate the ratio Z?_/Z)+. 
The Field-Feynman fragmentation functions used in the present analysis 
provide a good representation of the data from both EMC and HERMES 
(preliminary) (see Fig.Q) i.e. are close to those used in extraction of the d — u 
asymmetry by the HERMES collaboration. One should also remember that 
the effect of elastic p^ production is model dependent, being generally larger 
in the hard production mechanism than in the Regge model. If we take our 
Regge result at face value we argue that the d — u asymmetry extracted by 
the HERMES collaboration is rather overestimated. In the present analysis 
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data are from 1^. 
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we do not attempt to correct the HERMES data for the effects discussed 
here. This requires a separate analysis including efficiencies of the HERMES 
apparatus as well as knowledge of the many cuts used in their analysis. 

6 Conclusions 

Extraction of parton densities is one of the main goals of high-energy physics. 
It was proposed some time ago how to use semi-inclusive production of 
charged pions to determine both unpolarized and polarized parton densities 
in the nucleons. Recently this idea was put into practice in both cases. 
Such analyses assume implicitly the validity of the quark-parton model. In 
a recent work of two of us Q we have shown a breaking of the parton 
model in inclusive DIS at photon virtuality Q^ as large as 5-7 GeV'^ which is 
bigger than commonly perceived. The modern experiments analyzing semi- 
inclusive production of pions are performed in a similar range of Q^ . In 
general the situation in semi-inclusive reactions can be even more complex 
and subtle. In the present analysis we have made a first attempt to determine 
the nonleading mechanisms. 

We have estimated a few effects beyond the quark-parton model which 
may influence the extraction of the d — u asymmetry from semi-inclusive 
production of pions in DIS. 

Based on the analysis of hadronic data we have found that the inter- 
action of the resolved hadron-like photon with the nucleon may lead to an 
artificial enhancement of the measured d — u asymmetry in the region of 
small X. Next, we have investigated the elastic production of p^ mesons by 
a virtual photon on the proton and neutron targets based on two different 
models. Unequal cross sections for proton and neutron targets also lead to 
an artificial modification of the d — u asymmetry extracted based on QPM 
formulae. The effect found is opposite to the effect due to the resolved pho- 
ton component. These two effects cancel only in a narrow range of z. The 
net effect turned out to be z-dependent invalidating somewhat averaging in 
z as done recently in Q . 

We suggest that instead of averaging over a broad range of z one could 
try to select the region of z {x- and Q^-dependent) where the influence 
of nonpartonic effects is small. Unfortunately this can only be done at 
the expense of lowering the statistics considerably. An optimal choice of 
kinematical cuts in x, Q^ and z requires a more detailed study. Clearly, 
increasing of Q^ looks helpful. This could be realized by HERMES or by an 
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as yet unproposed experiment for COMPASS at CERN. 

Nuclear effects, even small ones, may also cloud the extraction of the 
true d — u asymmetry. To our knowledge there is no reliable estimate of 
such effects for semi- inclusive production of pions. 

Although in the light of the present analysis the precise direct extraction 
of (i — n is rather difficult, the semi-inclusive data can be used for tests of 
parton distributions, in particular the difference between d and -u, provided 
nonpartonic and nuclear effects are understood and included in the analysis. 

Finally we would like to point out that some of the effects discussed in 
the present paper may also influence the extraction of the polarized quark 
distributions from semi-inclusive production of pions in DIS. This will be a 
subject of a separate analysis. 

7 Appendices 

A Parametrizing the np -^ nX spectra in the pion hemi- 
sphere 

The mechanism of inclusive pion production in hadronic reactions is in gen- 
eral not well understood. It is believed that rather soft processes dominate. 
Some progress was made recently in studying inclusive pion production in 
polarized scattering (for a review see [^]). In proton-proton scattering the 
large-xj? region depends on the flavour structure of the outgoing pion. This 



part of the spectrum was explained e.g. in the recombination model |67], 
fusion model ^8| and recently in the meson cloud model [g^, ^, 0]. The 
central (mid-rapidity) part of the spectrum seems to be flavour independent 



The inclusive spectra of pions in pion-nucleon scattering are, in general, 
even less well understood. For beam-like pions the large-Xi? part of the 
spectra seems to be dominated by diffractive processes due to pomeron or 
reggeon (/ or p^) exchanges. For both beam-like and beam-unlike spectra 
one may expect a nonnegligible contribution from the decay of p mesons 
produced in peripheral processes vr-l-p— > p + X -^ (vr + vr) + X, dominated 
by pion exchange. 

The most complete experimental data for the vr^p -^ tt^X reactions were 
collected by the ABBCCHW collaboration at the CERN hydrogen bubble 
chamber [^6|. A detailed analysis of the vr+p — > ir^X and 7:~p — > n^X 
two-dimensional spectra p^ combined with a general understanding of the 
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reaction mechanism have shown that the spectra of all 4 reactions (4 x 10 
= 40 spectra) can be represented by the following six-component Ansatz: 

dxFdp'l V VmaxJ OXp 

+ Ctr, ■ haUxp) ■ e-^--Pi 

+ C'^^ ■ e± • e-^-"Pi 

+ Cip ■ xf{1 - xfT'^ ■ e-^^P'i 

+ C^- • /,(xi.) • e~^-Pi , (41) 

where the maximal rapidity r]max = Vmaxip])- 

Each of the components above corresponds to a distinct physical mech- 
anism, the first three to central processes and the last three to peripheral 



processes. By analogy with pp collisions |68] we have assumed one universal 
(flavour independent) soft component and allow for different normalization 
of flavour dependent components, called here hard due to their transverse 
momentum dependence. We have found Bgoft = 8.5 GeV~^ and Bhard = 
3.0 GeV~^), consistent with characteristic slopes of soft and hard processes. 
The parametrization of the soft component gives a multiplicity rising with 
the energy in the entrance channel. Some models in the literature predict 
a growth of the flavour asymmetric part of the inclusive cross section with 
energy, some predict that it should stay almost constant. Therefore in the 
present paper we have tried different functional forms for the phenomeno- 
logical flavour dependent part called here the hard component: 

{("]_ _ 2;2 \Phard 
^ rimax ' dxp 

In the present schematic parametrization we have taken ajp = 1 and 
aR = 0.5, i.e. we have neglected t-dependence of the Regge trajectories. By 
comparison with the experimental data we found B = 5.0 GeV~^. The extra 
factor xp in front of the diffractive components was added to extrapolate 
smoothly down to xp = 0. We have found that peripheral production 
of p mesons due to pion exchange and their subsequent decay constitutes 
a nonnegligible source of charged pions. Inspired by the pion exchange 
model and consistently with the data we have fixed the relations between 
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Figure 18: An example of the quality of the parametrization {41). 
experimental data points were scanned from Fig.l in ^ 



The 



normalization constants for such different processes 



C.1 



c: 
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The functional form of fp{xF) has been taken from a schematic model calcu- 
lation and one normalization parameter Cp was fitted to the two-dimensional 
spectra, which is possible because this mechanism dominates the beam- 
unlike spectra at xp > 0.7. 

The nature of the phenomenological very central, very soft {Been = 20- 
30 GeV~^) component is not clear. It was introduced only to to describe 
the data. It is most probably associated with pions from the decay of non- 
peripheral p mesons. We have found empirically approximately the same 
relation for normalization constants Crin as for Ca' . 



In Fig, 15 we present the quality of our fit for transverse momentum 



integrated Xi?-distributions for all four reactions vr^p -^ ir^X. The results 
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are shown in terms of the invariant single particle structure function 1 3t ] 

fM=- r^j^dpi. (43) 

71" Jo Pmax dxpdp'l 

B p meson decay functions 

In order to calculate the decay function in the most general case one needs 
to include off-diagonal elements of the density matrix [^2[. In the present 
paper we shall neglect the probably small off-diagonal effects. Then the 
decay function f{z) depends on the helicity of the parent p^ meson. In the 
general case of broad resonance it can be calculated as: 



/a(Vp) = j dmpp{m,) J dn[h{9, <f>)6{z{9, 4>) - z^/p)] , (44) 

where fx{0,(p) = \Yix{9 , (/})\'^ is the angular distribution of pions in the rest 
frame of p and p{mp) is the spectral density. The momentum fraction of a 
pion with respect to the parent p meson is: 

p|K)±MlV).l(l + eos^), (45) 

'^ nip 1 

which is independent of mp. The last relation is due to the smallness of 
the pion mass and must be corrected in the case of soft pions. In the 
approximation used in the present paper the two integrals in (|4^ factorize 
and one easily gets 




z) for A = ±1 
1)2 forA = 0, 



2 . , . (46) 



where above z is used instead of z^/p for brevity. 

For the exclusive reaction 7*A^ — > (PN at high energy one has z = z^/^ ~ 
z^ip. For semi-inclusive production through quark hadronization the decay 
function ( ^4[ ) must be convoluted with the fragmentation function into the 
p meson Dq^p. Below we shall consider these two distinct cases of p meson 
production. 

In inclusive unpolarized production of p mesons one may expect approx- 
imately an equal population of different helicities due to the complexity 
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Figure 19: Ry as a function of the photon virtuality. The experimental 



data are from [75, \7q] while the parametrizations are from ,^/. The long- 
dashed and the short-dashed lines correspond to the 2- and ^-parameter fits, 
respectively. 



of the poorly understood hadronization process. Then the effective decay 
function, which is averaged over p meson hehcities, becomes 



/(z)«-[2/i(z) + /o(z)]« const 



(47) 



In the most general case of exclusive p meson electroproduction the an- 
gular distribution of pions can be obtained according to the formalism of 
Schilling and Wolf |7^. For sufficiently large energy, where s-channel he- 
licity conservation takes place, averaging over azimuthal angle, the effective 
decay function can be approximated as 

fi{z) + eRv{Q\W)fo{z) 



^^""^'^ l + eRviQ^W) 

1/2 

where the polarization parameter e = e(Sg^,VK, Q 



(48) 



i-jt/ 



l-y+y2/2 



measures 
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the degree of longitudinal polarization of virtual photons. The effective 
decay function (48) depends on Q^ due to the empirically known strong 



Q^-dependence of Ry = a{'yLN -^ p^N)/a{'~^TN -^ p^N). Some smooth 
dependence of Ry on VK or a; is not excluded a priori. We shall take the 
model parametrization of RyiQ'^, W) from |74| (the 4-parameter fit) which, 



as shown in Fig.^, adequately describes the experimental data from [^, | 
We show there also recent HERMES experimental data [£^ which lie below 
the parametrization. Since the HERMES data were taken on the ^He target 
this may be partially caused by poorly understood nuclear effects. 
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